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A review of the main results of detailed flow analysis in highly central and semi-central heavy ion 
collisions at SIS energies is presented in the first part of this paper. The influence of the mass of 
the colliding nuclei and centrality on the collective expansion and the information on the equation 
of state of compressed and hot baryonic matter is discussed. The second part is dedicated to a 
similar type of analysis, based on the behaviour of the average transverse momentum as a function 
of mass of different hadrons, at the other extreme of energy range, where free baryonic fireballs are 
produced. Information on the partonic and hadronic expansion, temperature and degree of thermal 
equilibrium in p+p and Au+Au central collisions at 200 A-GeV is presented. 
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I. INTRODUCTION 

As far as the new states of matter are obtained in labo- 
ratory using heavy ion collisions, it is mandatory to have 
as much as possible under control the finite size and dy- 
namical aspects. One of these is the dynamical evolution 
of the transient piece of matter produced in heavy ion col- 
lisions. The present experimental evidence supports the 
existence of collective expansion caused by the pressure 
gradients built-up during the first phase of the collision. 
Understanding this expansion one could aim to pin down 
the properties of the initial phase of matter, before ex- 
pansion and at the same time, to extract information on 
in-medium effects, equation of state or phase transitions. 

Results of a detailed analysis of the expansion prop- 
erties of hot and compressed baryonic matter at SIS en- 
ergies in highly central and mid-central heavy ion colli- 
sions will be presented in Chapters I and II, respectively. 
Chapter III presents results of a similar analysis in terms 
of collective expansion and temperature at RHIC ener- 
gies, in the free baryonic sector. Preliminary results seem 
to support the possibility to disentangle between the con- 
tributions coming from partonic and hadronic phases and 
to access information on the degree of thermal equilibra- 
tion of these two stages of expanding matter formed in 
central collisions at 200 A-GeV. Conclusions will be pre- 
sented in Chapter IV. 



II. HIGHLY CENTRAL SYMMETRIC 
COLLISIONS 

For the highly central Au + Au collision at 150 A-MeV 
incident energy [l[ , the mean kinetic energy of fragments 
divided by their mass, was observed to present a depen- 
dence as a function of atomic charge which deviates from 
the one typical for a thermal plus Coulomb scenario. This 
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experimental trend agrees rather well with a thermal mo- 
tion superimposed onto a collective velocity field which 
in the non-relativistic limit can be written as: 

< Et? n >= \A ■ m - < pj low > +|"T" (1) 

\m - < (3f low > is the flow energy per nucleon and "T" 
stands for the effect of non-explicit treatment of Coulomb 
effects in the above expression. This leads systemati- 
cally to an overestimated value of the real temperature 
T [2J. While this trend was qualitativelyreproduced by 
the QMD microscopic transport model [3], the absolute 
value of mean kinetic energy per nucleon is systematically 
underestimated for all fragments. The experimental val- 
ues were reproduced rather well by a hybrid model based 
on hydrodynamical expansion coupled with a statistical 
fragmentation at the break-up moment [J] . This supports 
the model prediction that different species originate with 
different probabilities as a function of position within the 
fireball and time, heavier fragments being emitted later, 
at lower temperature and lower expansion. This was con- 
firmed in the meantime by small angle correlations stud- 
ies for pairs of nonidentical reaction species @] , [1] . A de- 
tailed analysis of experimental FOPI-Phase II data based 
on the equation above was performed for three symmet- 
ric systems (Au+Au, Xe + Csl, Ni + Ni), different in- 
cident energies and two regions of polar angles in the 
center of mass system for highly central collisions QjIH- 
While in the forward polar region, 25° < 9 cm < 45°, all 
three systems show the same flow at all measured en- 
ergies, along the transverse direction and mid-rapidity, 
i.e. 80° < 6 cm < 100°, the collective expansion is lower, 
its dependence on the incident energy is not that steep 
and at the same incident energy the collective expansion 
increases with the mass of the colliding system. Calcula- 
tions based on IQMD model at 250 A-MeV showed that 
the nucleons emitted at 90° suffer in the average 3.5 col- 
lisions relative to 1.7 of those emitted at smaller polar 
angles. This tells that at 90° the observed effect is most 
probably coming from an equilibrated fireball, while at 
forward angles the Corona effects and uncertainty in im- 
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FIG. 1: Collective energy, E co n, as a function of incident 
energies in highly central Au+Au collisions (1% total cross 
section). Open circles 25° < 9 cm < 45°, full circles 80° < 
dcm < 100°. IQMD predictions are represented by open and 
full squares for hard equation of state (HM) and open and full 
stars for soft equation of state (SM), for the corresponding 
polar regions. BNV prediction using a soft equation of state 
is represented by dotted dashed line 



pact parameter selection play an important role. Exper- 
imental slopes obtained for highly central Au+Au colli- 
sions, 1% total cross section, selected using the ratio of 
transverse and longitudinal energies of all detected and 
identified charged particles (E rat = J2i E±,i/ J2i m 
the cm. system) presented in Fig. [T] 9], support the mi- 
croscopic transport model estimates, IQMD Q and BNV 
[HI based on soft equation of state. 



III. MID-CENTRAL COLLISIONS 

What else could one learn at these incident energies 
going to mid-central collisions ? As far as concerns the 
centrality dependence of collective expansion one has to 
be aware that an average on azimuth (Fig. [2]) automat- 
ically includes contribution coming from shadowing ef- 
fects due to spectator matter while the values extracted 
at 90° (Fig. [3]) are strongly influenced by the aspect ratio 
of the fireball. With this in mind, one can see a system- 
atic increase of flow energy with the centrality (A part es- 
timated in a sharp cut-off geometrical model) the effect 
being more pronounced at the higher measured incident 
energies. A detailed analysis of the collective expansion 
azimuthal distributions at mid-rapidity, different impact 
parameters and incident energies for Au + Au and Xe 
+ Csl has been done by the FOPI Collaboration [ill ]. 
The qualitative agreement between the predictions of the 
model based on hydrodynamical expansion coupled with 
a statistical fragment formation at the break-up moment 
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FIG. 2: Azimuthaly averaged E co u as a function of A par t for 
Au+Au collisions at 90, 120, 150, 250 and 400 A-MeV and 
80° < (9cm < 100°. 



[4| , under the hypothesis that the expansion starts at the 
maximum overlap and perfect shadowing of the specta- 
tor matter, let us to conclude that different regions of 
azimuth capture different periods of central fireball ex- 
pansion 0],[ll|. If this is the case, then the amplitude 
of the azimuthal oscillation of E co u is recommended as a 
sensitive observable to the equation of state of hot and 
compressed baryonic matter produced at these energies. 
Transport model calculations based on BUU code [l2j 
using momentum dependent mean fields (m*/m=0.79), 
in-medium elastic cross sections (a = (Tq tanh(er^ ree /ero) 
with (To = P~ 2 ^ 3 ) and soft (K=210 MeV, gray zone) or 
stiff (K=380 MeV, dashed zone) EoS are compared in 
Fig. [4] and Fig. [5] with the experimental results. One 
should mention that the light fragments (up to A=3) are 
produced in a few-nucleon processes inverse to composite 
break-up, relative to the general coalescence recipe used 
by microscopic transport codes. 

The measured relative yields are nicely reproduced 
by this model, especially at higher incident energies 
[2j. For comparison with the experiment, the calcu- 
lated azimuthal distribution has been smeared according 
to the measured reaction-plane dispersion values. The 
azimuthal distribution of the average kinetic energy of 
different species and of the collective energy per nucleon 
extracted from these were fitted with: 
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FIG. 3: E co u as a function of A par t for Au+Au collisions at 
90, 120, 150, 250 and 400 A-MeV, 80° < 8 cm < 100° and 
72° < <f) < 108° ,252° < <j> < 288° (same symbols as in Fig. [2J 
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FIG. 4: E\ in and AEkin as a function of A par t, for Z—l 
(A=l,2,3) fragments, Au+Au at 400 A-MeV. The experimen- 
tal results are represented by triangles, while the BUU results 
are represented by gray zones for soft EoS and by dashed zones 
for stiff EoS, respectively. 



respectively. As one could see in these figures, E\ in for 
Z=l (A=l,2,3) fragments (Fig. gj) and E° oll (Fig.© are 
very little sensitive to the equation of state, both EoS pa- 
ramctcrizations showing quite good agreement with the 



data. As far as concerns AEkin and AE r 



the calcu- 



lations with the soft EoS reproduce the overall trends of 
the experiment while the calculations with the stiff EoS 
overestimate significantly the AEkin and AE co u values 
at higher and lower centralities, respectively. The results 
presented above support the conclusion that the equation 
of state of baryonic matter at densities of about 2p and 
temperatures of about 50-70 MeV is soft. This conclusion 
seems to be also supported by kaon production in heavy 
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FIG. 5: E® oll and AE co u as a function of A par t, for Au+Au 
at 400 A-MeV. The symbols have the same interpretation as 
in Fig. H 



ion collisions as it was shown by KAOS Collaboration 



.13]. 



IV. TOWARDS BARYONIC FREE MATTER 

The experimental information obtained at AGS, SPS 
and RHIC energies confirmed that the particle distribu- 
tions reflect mainly the conditions reached by the system 
in its final state. Therefore, information originating from 
earlier stages is mandatory in order to conclude whether 
the system passed through a partonic phase and extract 
information on such a phase. In central collisions, the 
final collective transverse, azimuthally isotropic flow, cu- 
mulates any collective flow generated during the evolu- 
tion of the fireball, i.e. from the partonic and hadronic 
phases, the yield ratios reflecting mainly the statistical 
nature of the hadronization process. This was nicely 
shown at RHIC, by the STAR Collaboration by an- 
alyzing the transverse momentum distribution for com- 
mon and rare particles using the blast wave model. The 
extracted kinetic freeze-out temperature and the flow ve- 
locity show a clear dependence as a function of central- 
ity while the chemical freeze-out temperature, extracted 
from particle yield ratios using statistical model stays 
constant, at about 170 MeV. At the same time the data 
seem to indicate a sequential freeze-out of particles, hy- 
perons decoupling from the system earlier at tempera- 
tures closer to the chemical temperature value. As we 
have seen in the first part of the present paper, a useful 
observable for a detailed study of collective flow would be 
the mean kinetic energy as a function of mass of different 
species. In the following we shall present the results of a 
similar analysis, this time in terms of the average trans- 
verse momentum (p t ) as a function of hadron mass for 
different colliding systems and energies. A careful exam- 
ination of the experimental (p t ) values [l5| as a function 
of mass for 7r ± , K ± , p and p evidences that the slope 
increases going from p+p to Cu+Cu and Au+Au at the 
same incident energy, i.e. 200 A-GeV. Although less pro- 
nounced, similar trend is observed as a function of inci- 
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dent energy for a given system [l6|]. In order to extract 
quantitative information we used the blast wave model 
for calculating the (p t ) as a function of mass: 



„ ^ Jo Ptf(pt)dp t 
< Pt >= 755 — 77-VT- 
Jo PtJ{Pt)dpt 



where: 
fiPt) 



rdrmt Iq 



Ptsinhp 
T 



nit cosh p 
T 



and 



p = tanh 1 (3 r 



A-(r)=A(£) 



(4) 



(5) 



(6) 



(7) 



The temperature T and expansion velocity (3 S were the 
free parameters used to fit the experimental (pt) as a 
function of mass [l7j]. Within this ansatz (3 = |/3 S . The 
results are presented in Fig. [5] The left part shows the 
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FIG. 6: Left part: the temperature T and expansion velocity 
/3 as a function of total mass of the colliding system at 200 
A-GeV; Right part: the temperature and (3 for Au + Au as a 
function of center of mass energy per nucleon. 

temperature and expansion velocity [3 as a function of 
total mass of the colliding system at 200 A-GeV and 
the right part the temperature and (3 for Au + Au as 
a function of center of mass energy per nucleon. The 
temperature drops by ~ 20 MeV and the expansion ve- 
locity is increasing from p+p towards heavy systems. Al- 
though the decrease in temperature is moderate it is con- 
sistent with the interpretation that stronger expansion 
observed in heavier combinations cools down the fire- 
ball and therefore the final, kinetic break-up tempera- 
ture is lower. Let's look now to {pt) as a function of 



mass for all particles measured by the STAR Collabo- 
ration at 200 A-GeV for p+p and Au+Au collisions [H[ 
presented in Fig. [7] The two lines represent the results of 
the fits using the above expressions for p+p ( dark line) 
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FIG. 7: Mid rapidity average transverse momentum as a func- 
tion of mass of different particles measured by STAR Collab- 
oration. 



and Au+Au (gray line). The corresponding temperatures 
and (3 values are: T=111.6±23.8 MeV and /3=0.39±0.06 
for p+p and T=I09.8±16.5 MeV and {3=0. 50+0.04 for 
Au+Au. While the freeze-out kinetic temperature seems 
to be similar for p+p and Au+Au, the expansion is 
much more violent in the Au+Au case. If for Au+Au 
one considers only tt^, K^, K* , K®, p, p, d, d parti- 
cles, the obtained temperature is T=98.7±19.5 MeV and 
/3=0.54±0.04. It is well known that for understanding the 
particle production in high energy and nuclear physics, 
many authors used Tsallis statistics [ijj]. Tsallis' gener- 
alization of Boltzmann-Gibbs extensive statistics, based 
on the definition of a q-deformed entropy functional, is 
supposed to be adequate for describing systems charac- 
terized by memory effects and long range-interactions. 
The Boltzmann-Gibbs statistics is recovered in the limit 
q— >1, therefore (q-1) is interpreted as a measure of the 
degree of non-equilibrium, the temperature T being in- 
terpreted as the average temperature. If we replace in 
Eq. (4), i(pt) with the distribution corresponding to the 
blast wave model in which the Tsallis statistics has been 
implemented 0|, [U: 
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and fit the same experimental data with the new expres- 
sion, the results presented in Fig. [5] are obtained. For 
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FIG. 8: Left part: the temperature T, expansion velocity /3 
and q parameters as a function of total mass of the colliding 
system at 200 A-GeV; Right part: the temperature, /3 and 
(q-1) parameters for Au + Au as a function of center of mass 
energy per nucleon. 



If all species are considered, as in Fig. [7J the results of 
the fit are summarized in Table U where the blast wave 
models with Boltzmann-Gibbs (BGBW) and respectively 
Tsallis (TBW) statistics are compared. (3 is almost zero, 
therefore no expansion is built up in p+p collision at 200 
A- GeV and the degree of non-equilibrium is about a fac- 
tor of four larger in p+p relative to Au+Au. For Au+Au 
the extracted temperature is about 20 MeV lower and (3 
is within the error bars similar with the one obtained 
using the blast wave model with Boltzmann statistics 
scenario. If one considers separately common particles, 
i.e. hadrons with larger interaction cross section, the 
results are presented in Table [TTJ The temperature de- 



TABLE II: 





Au + Au 


Au + Au 




BGBW 


TBW 


T [MeV] 


98.7±19.5 


79.05±0.04 


P 


0.54±0.04 


0.53±0.0005 


q 


1.0 


1.0175±0.0018 



creases, (3 slightly increases and the deviation from Boltz- 
mann statistics is reduced. If we consider only the hy- 
perons and J/'!', as it is shown in Table ITTTI a tempera- 
ture around 200 MeV, j3 ~0.3 and larger deviation from 
global equilibrium is obtained. This suggests that strange 



TABLE III: 



TABLE I: 



System 


p + p 


P + P 


Au + Au 


Au + Au 


Model 


BGBW 


TBW 


BGBW 


TBW 


T 


111.6+23.8 


78.86+10.13 


109.8+16.5 


86.8 ±1.54 


P 


0.39+0.06 


0.027+0.10 


0.50+0.04 


0.48+0.04 


q 


1.0 


1.087+0.002 


1.0 


1.025+0.003 



p+p the boosted Tsallis scenario gives a negligible ex- 
pansion velocity. As far as concerns the q parameter, 
one could observe a strong decrease as a function of total 
mass of the colliding system, large deviation from Boltz- 
mann statistics, (q-l)~0.085, for p+p to relative small 
value, 0.005, for Au+Au being observed. The value for 
Au+Au stays almost constant, ~0.005 as a function of 
incident energy. One should mention that a similar trend 
as a function of incident energy but for larger (q-1) values 
was observed in e + e~ collision, obtained from analyzing 
the p t spectra in terms of simple, not boosted, Tsallis 
distribution 12211. 





Au + Au 




TBW 


T [MeV] 


198.0+7.6 


13 


0.32+0.012 


q 


1.0247+0.0043 



and heavy flavour hadrons keep the characteristics of ex- 
pansion at the hadronization moment, characterized by 
lower values of /3, higher temperature and based on Tsal- 
lis statistics interpretation, not fully equilibrated. After 
hadronization, strongly interacting particles, within still 
highly dense environment continue to build up expan- 
sion, cooling down the system and approaching a global 
equilibrium. 

The Minuit package was used to perform the least \ 2 
fit. The quality of the fit was reasonably good as can be 
seen in Fig. [5J where 1/x 2 is represented as a function of 
two pairs of parameters for the result of the fit listed in 
Table III. 
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FIG. 9: a) 1/x 2 as a function of temperature and /3; b) 1/x 2 
as a function of temperature and q for the result of the fit 
listed in Table III. 

V. CONCLUSIONS 

In this paper we tried to present a short review of 
the main results of detailed flow analysis at incident en- 
ergies where compressed and hot baryonic fireballs are 
produced. The experimental data support model predic- 
tions that different species originate with different prob- 
abilities as a function of position within the fireball and 
time, heavier fragments being emitted later, at lower tem- 
perature and lower expansion. The best agreement be- 
tween experiment and theoretical estimates based on mi- 
croscopic transport models is obtained if a soft equation 
of state is considered. 

Similar type of analysis extended at the other ex- 
treme of incident energies, where free baryonic matter is 



produced, shows that for central collisions, azimuthally 
isotropic flow of different species could be used to pin 
down the relative contributions coming from partonic 
and hadronic levels. The kinetic temperature and expan- 
sion velocity extracted from experimental < p t > depen- 
dence on the hadron mass using boosted Boltzmann and 
Tsallis expressions support the conclusion that strange 
and heavy flavour hadrons carry the information from 
the partonic stage characterized by higher temperature 
and lower expansion. The results based on boosted Tsal- 
lis distribution (within the limits of its applicability and 
interpretation) indicate high degree of non-equilibrium 
and missing expansion in pp collisions at RHIC energy. 
For the Au+Au collision similar analysis for strange and 
heavy flavour hadrons gives a higher temperature, less 
violent expansion and higher degree of non-equilibrium 
of the initial, deconfined matter produced at the highest 
RHIC energy. Common hadrons, strongly interacting af- 
ter hadronization, continue to build up expansion, cool 
down the system and bring it towards a global equilib- 
rium. 

It is a real challenge for the near future experiments 
at LHC to use the potentiality of flow phenomena to ex- 
tract new physics in an energy domain of about 30 times 
larger than the one attainable at RHIC. Preliminary re- 
sults based on Monte Carlo simulations show that the 
ALICE experiment is able to reconstruct with high accu- 
racy the p t distribution for most of the particles making 
the type of analysis presented in this paper feasible and 
promising. 
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